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p I I Abstract 

T^ , In this paper, we derive a non-linear version of the Feynman-Kac formula for the 

C^ I solutions of the vorticity equation in dimension 2 with space periodic boundary condi- 

tions. We prove the existence (global in time) and uniqueness for a stochastic terminal 
value problem associated with the vorticity equation in dimension 2. 

> 

2i : 1 Introduction 

m 

The Feynman-Kac formula, in its original form derived from the idea of path integration in 
'nI" . Feynman's PhD thesis (which is now available in a new print [7]), is a representation formula 

for solutions of Schoingder's equations, and in the hand of Kac, is an explicit formula written 
in terms of functional integrals with respect to the Wiener measure, the law of Brownian 
motion. 

K^ • Bismut [3], Pardoux-Peng [TT] and Peng [12], by utilizing Ito's lemma together with Ito's 

\-{ . martingale representation, have obtained an interesting non-linear version of Feynman-Kac's 

■ - -' formula for solutions of semi-linear parabolic equations in terms of backward stochastic differ- 

ential equations (BSDE). The goal of the present paper is to derive a Feynman-Kac formula 
for solutions of the Navier-Stokes equations in the same spirit of Bismut and Pardoux-Peng 
[TT] , and to study the random terminal problem of the stochastic differential equations as- 
sociated with the vorticity equations. 

The main idea contained in [3], [TT] may be described as the following. Let u{t,x) = 
{u^(t,x), ■■■, u"^(t,x)) be a smooth solution to the Cauchy initial value problem of the 
following system of semi-linear parabolic equations 

d 

—u'-uAu' + f{u,Vu) = 0,u{0,x)=uo{x) in R'^ (1.1) 

where i = 1, ■ ■ ■ ,m, and z/ > a constant. Let B = {B^, ■ ■ ■ , B'^) be the standard Brownian 
motion on a complete probability space {Q, J-", P), x G M'^ and T > 0. Let us read the solution 



u along Brownian motion B. More explicitly, let Yt = u{T — t, \/2vBt + x) for t G [0, T] and 
Zt = VuiT — t, \/2uBt + x), and apply Ito's formula to u and B to obtain 



YT-Yt = j^ f{Y,,Zs)ds + V2^ /f Z, ■ dB, 
Yt = Uo{Bt). 



1.2) 



In literature, (11. ip may be written in differential form 

dY = f{Y,Z)dt + y/2l^Z ■dB,YT = ^, (1.3) 

where the arguments s, t etc. are suppressed if no confusion may arise. The differential 
equation above is an example of backward stochastic differential equations, where the termi- 
nal value It = ^ is given. The function / appearing on the right hand side of (11.31) is called 
the (non-linear) driver. 

Pardoux-Peng ^^ made an important observation. If the non-linear driver / in BSDE 
(II. 3p is globally Lipschitz continuous, then there is a unique adapted solution pair (Y, Z) 
satisfying (II. 3p for a random terminal value ^ G L^(r2, J-2-,P), which is not necessary in the 
form of uq{Bt)- The solution u and its gradient Vu in turn can be represented in terms of 
{Y,Z). This representation may be considered as a nonlinear extension of Feynman-Kac's 
formula to semi-linear parabolic equations. 

More recently, Kobylanski [8], Delarue [6], Briand-Hu |1], Tevzadze [13], and etc. have 
extended Pardoux- Peng's result to some BSDEs with non-linear drivers of quadratic growth. 
These papers however mainly deal with scalar BSDEs only, which corresponds to semi- 
linear scalar parabolic equations. It remains largely an open problem whether the BSDE 
approach may be applied to non-parabolic type of partial differential equations. We study 
in the present paper a class of backward stochastic differential equations which arise from 
the vorticity formulation of the Navier-Stokes equations, hence provide Feynman-Kac type 
formula for solutions of the Navier-Stokes equations. 

Relations between the Navier-Stokes equation and foward-backward stochastic differential 
equations formulated in the group of diffeomorphisms were introduced in ^ . 

2 The vorticity equation 

Let us describe a class of (infinite dimensional) backward stochastic differential equations 
associated with the study of Navier-Stokes equations. 

The 2D Navier-Stokes equations (without external force) are the partial differential equa- 
tions which describe the motion of fluids 

g - uAu + u ■ Vu + Vp = 0, , . 

where u = (m^,u^) is the velocity field, u the viscosity constant and p the pressure. The 
mathematical study of the Navier-Stokes is interesting by its own, and even the simplest 
situation where the space periodic condition is supplied is of interest. 



Suppose that u{0,x) = ip{x) is a smooth vector field with period one, that is, ip{x + ei) = 
if{x) for all X G M^, where Ci = (1,0) and 62 = (0,1) the standard basis in M^. Then, 
the unique solution {u,p) to the 2D Navier-Stokes equation is smooth on (0, cxd) x M^ and 
periodic in space variables, so that u{t,x + ei) = u{t, x) and p(t,x + ei) = p(t, x) for all t > 0, 
a; eM^ z = 1,2. 

Let 

0;=^-^ (2.2) 

dx^ dx"^ 

be the vorticity of u, which is a scalar function in dimension two, and thus the evolution 

equation for a; is a scalar partial differential equation 

-^ - z/Aw + M ■ Vw = 0. (2.3) 

ot 

Equation (12. 3p is called the vorticity equation is equivalent to the Navier-Stokes equation. 
The relationship between the scalar function u and the associated vector field u is determined 
by the Poisson equations 

A.' = -^ aud ^u' = ^. (2.4) 

By (12.21) the average of uj{t, x) 

uj{t, x)dx = for all t > 0, (2.5) 

[0,1)2 

so that (12. 4p has a unique periodic (with period one) solution u = (u^,m^). We define linear 
operators Ki : u —^ u^ (where i = 1,2) and K : u ^ u hy solving the Poisson equations 
(12. 4p . where w is a real function with period one and mean zero. 

Let T^ = M? /7? be the 2D torus equipped with the standard metric and the Lebesgue 
measure. We may identify tensor fields in M^ with period one canonically with the corre- 
sponding tensor fields on T^. For example 

L\T^) = {fe LlM') ■■ /(• + e.) = /(•) for z = 1, 2} n ^^([0, 1)^). 

If / G L2(T2) then 

f{x) = J2 e'-^'^'^-^fik) (2.6) 

where 

/>) = f e-'''^'^''y^f{y)dy, k G Z' (2.7) 

"'[0,1)2 

is the Fourier transform of /, (■, ■) denotes the scalar product in Euclidean spaces. 
Lemma 2.1 (Green's formula) Consider the Poisson equation 

Ag = -f m T\ [ g{y)dy = 0, (2.8) 

./t2 



where Jj2 f{y)dy = and f G L^(T^). Then the unique solution of the problem Ii2.8\) is given 
by 

9i-)= E ^d^^^»- (2.9) 

Our first goal is to derive a probabilistic representation for u in terms of Brownian motion. 
To this end we set up the probability setting with which we are going to work with. Let 
B = [B^^B^) be a standard Brownian motion on a complete probability space (fi,^-', P), 
and define 

Y{w,t,x) = u{T-t,x + y/2^Btiw)), 

Z\w,t,x) = §,{T-t,x + V2^Bt{w)), 

duj 



Z\w,t,x) = §,{T-t,x + V2vBt{w)) 

for {w,t,x) E Q X [0,00) X M^. We will often suppress the random element w from our 
notations, and write Y{t,x), Yt or simply by Y for Y{w,t,x), if no confusion is possible. 

Let ijj = gfr — g^ be the vorticity of the initial velocity ip = Mq, and ^(x) = 'ijj{x+y/2vBT)- 
Then, it is clear that ^ is smooth and periodic in x (with again period one). According to 
Ito's formula 



^{x)-Y{t,x) = V2^ {Vuj{T-s,x + V2^Bs),dBs) 

(—^ + uAuj\{T-s,x + V2^B,)ds. (2.10) 

Now, by utilizing the vorticity equation (12. 3^ : substitute — ^ + z^Ao; by u ■ Vw to obtain 

nT nT 

i{x)-Y(t,x)= j {Z{s,x),X{s,x))ds + V2^ {Z{s,x),dB,) (2.11) 

where for simplicity we have set 

X(t, x) = u{T-t,x + V2vBt) 

which is continuous in t, smooth in x, and periodic in x. Next, we wish to rewrite X{t,x) 
in terms of Y and Z. To this end, it is a good idea to introduce some notions in Fourier 
analysis, and establish several notations which will be used in what follows. 

Let us apply Green's formula to the vorticity u of u. According to (12. 4p and (12. 9p . we 
have 

^ E ^e-v-T(^-)^(A:) (2.12) 



u^(x) 



2 ^ , \k[^ 



and 



u^(x) 



J2 Ae2-v^('=.-)S(A;). (2.13) 



2 ^^ |A;|2 

fc=(fci,fc2)ez2,fc^o 



In other words 



^ / 1 u ^ I 1 /j. 



u{k), fc ^ 0. (2.14) 



Hence 



X^(t,x) = u\T-t,x + V2^Bt) 



^ E 



2 ^^ |A;|2 



^ 27rv'^(A:,x+v^St> 



a;(T-t, Ol^). (2.15) 



On the other hand 



Yit,-){k) = I e 

'[0,1)2 



^2nV^(k,y) 



io{T-t,y + V2vBt)dy 



= e^-^^'^'^^^^ / e-2-v^<'='S^)a;(r - t, y)dy 

i[0,l)2+v^Bt 

J[0,l)2 

= e2'^^^<'='^^*>a;(f^,-)(A;) 

here the third equahty follows from the fact that u is periodic. Substituting the above 
equality into (12.151) to obtain 



x\t,x) = y^ J2 



k 



2 ^^ |A;|2 



^ 2n^/^{k,x) 



Yit,-m, 



and 



X^{t,x) 



^ E 



2 ^^ |A;|2 



^1 e2-v^('=.-)F(^)(A;). 



(2.16) 



(2.17) 



By our definition of linear operators Ki, the previous equations (I2.16[ I2.17P may be written 
as 

X^{t,x)=Kj{Y{t,-)){x) Va;GR^J = l,2. (2.18) 

Thanks to (I2.18P we may finish our computation for Y as following. According to (12. lip 

i{x)-Y(t,x) = V^ [ {Z{s,x),dBs)+ [ {Z{s,x),K{Y{s,-))ix))ds (2.19) 



where x runs through 



3 Feynman-Kac formula for the vorticity 

The preceding stochastic integral equation f l2.19p may be put in a differential form 

dY = {Z,K{Y))dt + V2^{Z,dB), Yt = i (3.1) 

where the time space variable x, for simplicity, is suppressed. The initial value problem to the 
vorticity equation f l2.3p is transferred to a terminal value problem to the stochastic differential 
equation ( 13. ip within the formulation of BSDEs. In order to derive a non-linear version of 
the Feynman-Kac formula for the vorticity w, we need to study the infinite dimensional 
BSDE (I3A|) . 

BSDE (13.1 p possesses two features which make it difficulty to study. First, the stochastic 
equation (13. ip must be solved in a function space, so it is an infinite dimensional stochastic 
differential equation (with finite dimensional noise). Second, the non-linear term in this 
BSDE is quadratic in Y and Z, which is the origin of all difficulties. There are few results 
in literature about this kind of backward stochastic differential equations. 

Let B = {B^,B^) be a standard Brownian motion on a complete probability space 
{fi,J^,F). Let J-'° = a{Bs: s < t} and {J-'t)t>o be the completed continuous filtration 
associated with (J^°)t>o. Let O and V be the optional and predictable cr-fields on fi x [0, oo), 
respectively. Let Q = O x i3(M^) and V = V x B{R^) be the optional and predictable a- 
algebras on f2x [0, oo) xM^. A Q- measurable (resp. P-measurable) function on f2x [0, oo) xR^ 
is called a optional (resp. predictable) function, or called an R^-valued optional (resp. 
predictable). We may similarly define O x i3(T^) and V x i3(T^) which are identified 
with elements in the O x i3(]R^) and V x i3(]R^) respectively which are periodic in the space 
variables with period one. 

In order to derive a non-linear version of the Feynman-Kac formula for the vorticity 
u we need to prove the existence and the uniqueness of solutions to (13. ip subject to the 
given terminal value ^. Actually we will do this for a general terminal value ^ which is not 
necessary in a form of iplB^). 

We will assume that ^ is a bounded random function on ^2 x T^ which is J-'t ® i3(T^) 
measurable, and furthermore, we assume that for every w & Q, ^{w,-) G ly^'^(T^), and 
Jj.2^{w,y)dy = for all w E fl. In particular, according to the Sobolev imbedding, x — )■ 
$,{w,x) is continuous on T^ for every w G Q. 

By a solution to BSDE (13. ID we mean a pair of P-measurable stochastic processes Y and 
Z, such that: 

1) For all X G T^, {w,t) — )■ Y{w,t,x) is continuous semimartingale, and for all {w,t) G 
fix [0,T],F(M;,t,-) gL2(T2). 



2) For every a; G T 



T 

2, 



E / \Z{t,x)\^dt <+oo 







so that the Ito's integral J^^{Z{-,x),dB) is a square integrable martingale for every x. 



3) It holds that 

Y{t,x)=^ix)- f {Z{s,x),K{Y{s,-))ix))ds-V2^ f {Z{s,x),dB,) 

almost surely on f2 x T^, for t G [0, T]. 

Now we are in a position to state our main result. 

Theorem 3.1 Under above assumptions on the terminal value ^, there is a unique solution 
pair (Y, Z) to BSDE (EJP such that 

1) Y is hounded onVl x [0,T] x T^, and 

2) For almost all x G T^, the ltd integral JQ{Z{-,x),dB) is a BMO martingale, and 



T 

2, 



Tt\ < +00 



ess sup E < / 1 1 Zs 1 1 (is 

[0,T]xQ \jt 

where || ■ || denotes the l?-norm on T^. 

In particular, by applying Theorem 13.11 to ^ = '^[Bt) where yj = V x uq is bounded, 
C^ on T^, then Y{t,x) = uj{T — t, ^j2vBt + x), where (F, Z) is the unique solution pair of 
(13. ip with terminal ^ = ^{Bt) and uj is the solution to the vorticity equation (12. 3 p subject 
to the initial value u(0, ■) = Uq. Y may be regarded as the probabilistic representation for 
the vorticity uj. 

The proof of Theorem 13.11 relies on two important technical facts. The first is the 1?- 
estimate for the linear operator K^ and the secomd is a maximal principle formulated in 
terms of backward stochastic differential equations. 

4 Several technical estimates 



In order to prove the main result Theorem 13. H we need several a priori estimates. 

4.1 Apriori estimates for K 

Let us recall the definition of K. Note that we identify tensor fields in the torus T^ with the 
tensor fields on M^ with period one along each space variable. For k G Z_|_ and q > 1 the 
Sobolev space 

W'^'^'iT^) = {/ : d"f G Ll^iM.^) n ^"([0, 1)^) for |a| < k 
and/(- + eO = /(-) for ^ = 1,2} 

together with the Sobolev norm 

1/9 

aeZ2,|a|<fc 



where 1 1 ■ 1 1 ^ is the L'^-norm over T^ , that is 

l/r = / \fix)\'dx 

T2 / \^ [0,1)2 

If g = 2 then we use 1 1 ■ 1 1 instead of 1 1 ■ 1 1 2 for simphcity. 

According to Sobolev's embedding theorem, ly^'^(T^) M- C"(T^) for some a G (0, 1), so 
any element in ly^'^(T^) has a unique continuous representation. 

If / G L^(T^) such that /fq ]^^2 / = 0, then Kj{f) = Qj are the unique solutions (with 
period one) such that L ^-.2 Qj = 0, solving the Poisson equations 

Ag^ = -|^, Ag, = |^ on T^. (4.1) 

0x2 OXi 



By definition, if a = (ai, 0:2) G Z+ x Z+, then d"Kj{f) = Kj{d°'f), where d" stands for the 
partial derivative 4° ^2 for simplicity as long as d°'f G L^(T^). 
On the other hand 

df 



/ iV^jf = - / gjAgj = / ^1-— or - / ^2 

Jt2 jt2 Jt2 ca^2 Jt2 



9xi 



according to j = 1 or j = 2. Integrating by parts together with Cauchy-Schwartz inequality 
to the last integrals we deduce that 



Jt2 Y Jt2 y Jt2 

which yields that 

||Vi^,(/)||< 11/11, J = 1,2. (4.2) 

Let Ai > be the spectral gap for the torus T^. Since J^2-^j(/) = 0, according to the 
Poincare inequality 

ll^.(/)ll<^l|Vi^,(/)||<^||/||. (4.3) 



Therefore we have the following elliptic estimate (for more details see for example [T], [2] 
and [lOj). 

Lemma 4.1 There is a universal constant Cq > such that 

\\K,{f)\\,,2<Co\\m-l,2 

for any f G W^-^^^{T^) with /^^ / = 0, where A; G N. 

In particular, if / G ly^'^(T^), K{f) is a-Holder continuous. 



4.2 A maximum principle 

Let us formulate a probabilistic version of the maximum principle in terms of BSDE. 

Lemma 4.2 Suppose y is a continuous semimartingale such that 



Vt 



Vt- f {h,z)ds- f {z,dB) fortE[0,T], 



where ht is a hounded J^T-i^^o,surahle random variable, both z and h are W^ -valued predictable 
processes such that 

E / \z\^dt < cx) 
Jo 

and suppose that 

Rt = exp - / {h,dB) - - / \h\'^ds 
I Jo 2 Jo 

is a martingale up to T. Then \yt\oo < IvtIoo for all t G [0, T] almost surely. 

Proof. Define a probability Q with density R. Then P is equivalent to Q on J-^, and 
according to Girsanov's theorem Bt = Bt + J^ hgds is a standard Brownian motion, and 

yt-yT = - / {z,dB). 
Therefore yt = ¥P {yrlJ^t} so that \yt\oo < IvtIoo almost surely. ■ 

4.3 A linear BSDE 

Let us consider the following linear BSDE 

dY ={Z,h)dt+{Z,dB), Yt = ^ (4.4) 

with h E Q is a given optional process (valued in T^) such that for each {w,t) E ft x [0,T], 
h{w,t,-) gC(T2) ^^^ 

E \h{t,x)fdt<oo VxgTI (4.5) 

Jo 

^ is the terminal value: 

^EL°°iQx T^) n L~(fi, J^T, W^'\T^)). 

The linear equation (14. 4p may be solved for every x eT^, and in fact we may solve the linear 
BSDE 

dY{t,x) = {Z{t,x),h{t,x))dt+{Z{t,x),dBt), , . 

Y{T,x) = ax), ^ ' 

9 



by means of changing probability. More precisely, for each x G T^, since (14. 5 p holds, we can 
define a probability Q"^ on T^ by ^^ = R{T,x), where 



R(t, x) = exp 



- f {h{s,x),dBs)-]- I \h{s,x)\''ds 
. Jo "^ Jo 



If {Y{-, x), Z{-, x)) is the unique solution of (14. 6p . then, according to the Girsanov theorem, 
Y{-,x) must be a martingale under the new probability Q^', we therefore have 



F(t,a;)=E^^{e(x)|J-,} 



which implies that 



^■('-'-^^f!ly«w 



J't 



for {t,x) G [0,T] X T^. Therefore we have established the following 

Lemma 4.3 Suppose that ^ is W'^''^{T'^) -valued J^T-'measurahle random variable, and sup- 
pose that h is a W'^''^(T'^) -valued adapted stochastic process satisfying lli4-5\ ), then the unique 
solution to (R3P ^'^ given by 



F(t,x)=E{e(x)e- 



- f^{h{s,x),dBs)-^ f^ \h{s,x)\'^ds 



Tt 



(4.7) 



forte [0,T] X T2. 
It is clear that 



d,Y{t,x) = eU.^x) 



R{T, x) 



R{t,x) 



Tt 



+E<;^(a:)(-/ {d,his,x),dB,)- J^ (^h{s,x),d^h{s,x))ds\^^ 



Tt 



so we have the following simple fact. 



Corollary 4.4 1) If in addition ^ and h are bounded, then the solution Y is continuous in 
{t,x). 2) If in addition ^ and h have bounded derivatives in x, then so is Y. 



5 Proof of Theorem 13.1 



This section is devoted to the proof of Theorem 13.11 We will use the following convention in 
our proof. The elliptic estimates show that if / G W'''^{T^) then K{f) G W''+^'^{T^), thus 
if fc > 1, then, according to the Sobolev imbedding, K{f) has a Holder continuous version. 
Therefore, if / G W^''^{T'^) for k > 1, K{f) is always chosen to be its continuous version. 



10 



Let "H denote the collection of all bounded P-measurable stochastic processes Y on 
Q X [0, T] X T^ satisfying the following conditions: 

1) For each x G T^, Y{-, x) is a continuous seminiartingale (up to time T) on [Q, J-', Tf, P) 
whose martingale part M (with Mq = 0) is a BMO martingale, and Yt = C,- Moreover, for 
every {w,t) eilx [0,T], Y{w,t,-) e W^'\T^). 

2) Let the Ito representation of the martingale part 

M{t,x) = I {Z{t,x),dBt) 
Jo 

where Z is P-measurable. Then 

Tt )■ < +00. 



ess sup E < / 1 1 ^s 1 1 "^ds 
[o,T]xn {Jt 

Let y G "H, we define C{Y) = Y hj solving the following linear backward stochastic 
differential equation 



dY{t,x) = {Z{t,x),K{Y{t,-)){x))dt + V2v{Z{t,x),dBt),Y{T,x) = ^{x). (5.1) 
for every x G T^. Then F G "H. 

5.1 Apriori estimate for the density process Z 

If F G "H, in particular F is a bounded function on 1] x [0,T] x T^. ||V||oo denotes the 
essential bound of F on fi x [0, T] x T^. 

Suppose y G "H such that ||l^||oo < Ci- Define Y = C{Y), and Z is the density process 
of the martingale part of Y, that is, {Y, Z) by solving the following linear BSDE 

dY = {Z,K{Y))ds + V2^{Z,dB),YT = i (5.2) 

where \^{w,t,x)\ < Ci. By the maximal principle, \Y{w,t,x)\ < Ci. 

For simplicity, we will use E-^* to denote the conditional expectation E{-|J-'(}. 
By Ito's calculus, 

fT rT 



\Yt\^ = |^|2_2z// \Z\^ds-2 I Y{Z,KiY))ds 
-2v^ / Y{Z,dB) . 
First take conditional expectations, to obtain that 



11 



\Yt\^ + 2v¥.^' / \Z\^ds = E-^*|^p-2E-^* / Y{Z,K(Y))ds 
Jt Jt 

cT 



< Cf + 2CiE^' f \{Z,K(Y))\ds 



then integrating over T^ and using the estimate from the maximum principle, we have 



iFtlr + 2z/E- 



J^t 



< Cl + 2CiE 



Tt 



ZW^ds 
\{Z,K{Y))\ds 



t JT2 
T 



< C't + 2CiE-'^ I \\K{Y)\\\\Z\\ds 

2 , r^^Tt I e\\K{Y)\\^ + -\\Z\ 



< Cf + CiE- 
Jt 

< Cl + eCiCoE^' 
for every e > 0. Recall that 



ds 



T 



iFirrfs 



e 



\Z\\^ds 



\Z\\qj^q = ess sup E * / \\Z\\ ds . 

nx[Q,T] Jt 



It follows that 



I7IP < *^^ 



1 
e 



Choose £ = — , we obtain 



Ci + TeC^ oq + - 1 1 ^1 1 BMO 



Ci 



(5.3) 



V 



That is, the norms ||i^||oo and ||^||_bmo are uniformly bounded, depending only on z/, Ci, 
Co and T. 

5.2 Contraction property 

Let q; be a real number to be chosen later, and consider Y^ = e°'^Yt and Y^ = e°'^Yt. Then, 
according to integration by parts 

rff" = (Z, K{Y''))ds + V2^(Z°, dB) + aY'^dt. 
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Denoting SY" = F" - Y"^ and 5Z'^ = Z" - Z"^. Then 

where 

$, = (Z„ir(F,"))-(Z;,ir(Fr)). 

It follows that 



|(5Fi"|' = -2z/ / \5Z''\'ds-2a \5Yyds 

-2 I {5Y'')Ms - 2V2^ f {SY''){SZ'',dB) 
by taking conditional expectation given the information up to J-t to obtain 

-2E-^* /" ((5F")<l>rfs. 
Now integrating over T^, to obtain 

||5i;"||2 = -2z/E-^* /" ||5Z"|pds-2aE-^* /" ||5f"|prfs 



-2E-^* /" /" ((5r°)$rfs (5.4) 

Let us write for simplicity 

J(t) = ||5F°||2 + 2z/E-^* /" ||(5Z"||2ds + 2aE-^* /" WdY'^W^ds. 
Then (15 ■4p implies that 

J{t) = -2E-^* /" /" ((5f°)<l>cis 

< 2E-^* f ||5F"||||$||rfs 

E-^* /" WSY^W^dsY (e^' f IM^dsY (5.5) 

which yields that 

< 2aE-^' / ||(5y"|Prfs + — E-^' / ll^lPrfs. (5.6) 

Jt 2a Ji 



13 



Let us now consider the last integral appearing on the right-hand side of (15. 6p . It is clear 
that 



T) 



ii$.ii = \\z,-K{Yn-z:-K{Y:' 

= \\Z,.K{6Yn + 6Z:.K{Y:)\\ 

< \\K{6Yn\\\\Z.\\ + \\K{Y:)\\\\SZ:\\ 

< Co\\6Y:\\\\Z,\\ + Co\\Y:\\\\6Z:\\ 

plugging into (15.61) we conclude that 

J{t) < 2aE^^ r \\6Y'^\\'ds + ^°^' '^ + ^f °^' \\6Y'^\\l 
Jt « t^ 



+r£nE-^* / Wbz^W'ds (5.7) 

a Jt 

where we have used the uniform bounds 



\y\U < Ci and \\Z\\bmo < —\/iy + TCoCf . 



Choose a > such that 



then (15. 7p yields that 



ClClv^TC^Cl ^ 1 ClCl ^ V 
a v^ ~ 16' a ~ 4 



||(5F"|U + ||5Z"||iJAfo<^||5F^ 



Theorem 5.1 There zs a > such that, C is a contraction on Ti under the norm 

1 1 VI I 1 1 V^l I _J_ 1 1 7" 1 1 

N-*^ lla.oo — II-' ||oo + ||^ WBMO 

where Zf = e'^^Zf and Z is the density process of the martingale part ofY. 

We are now in a position to complete the proof of Theorem 13. II The sequence of Picard's 
iteration is constructed as the following. Begin with 

Yo{t,x)=E{^{x)\J^t} 

(here we mean the continuous version of the optional projection of ^) and Zq is the den- 
sity process of Yq with respect to the Brownian motion determined by Ito's martingale 
representation. Since ^ e ly^'^(T^), so Yo{t,-) G ly^'^(T^) for all t almost surely. Define 
Yn+i = C{Yn) for ra = 0, 1, 2, ■ ■ ■ . Then Lemma [4.31 implies that all Y^ &% and in particular 
(t, x) — )■ y„(-,t, x) are continuous almost surely, so that 

f{\Yn{t,x)\ < Ci for all (t,x,n) G [0, T] x T^ x N} = 1. (5.8) 

Theorem 15. II implies that {^n} is a Cauchy sequence under the norm || ■ \\a,oa for some a > 0, 
and therefore has a limit Y which is a solution to (13. ip . 
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